The anomalous metallic ferromagnetic state of Sr doped manganites. 
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We deduce a model relevant for the anomalous metallic state of Sr doped manganites at low 
temperatures within the ferromagnetic phase. It provides a natural explanation to several anomalous 
features observed experimentally, such as the vanishing Drude contribution in optical conductivity, 
the pseudo-gap in the density of states, an the unusual dispersion observed in photoemission. 
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The fascinating physics of the perovskite manganite 
systems is continuously attracting the interest of ex- 
perimentalists and theorists alike. Initially, most of 
the attraction focused on the colossal magnetoresistance 
(CMR) effect, due to its potential technological appli- 
cations. Nevertheless, these systems brought us other 
remarkable and perhaps more puzzling surprises as, for 
instance, the strange nature of the metallic state at low 
temperatures. Optical conductivity experiments in Sr 
doped manganites La x Sri_ x Mn.03 by Tokura et al. |Q, 
have demonstrated the presence of an extremely small 
Drude part with a strong mid-infrared contribution dom- 
inating the low frequency response. The dc-resistivity, 
on the other hand, shows conventional T 2 metallic be- 
haviour, but its coefficient A has a strong doping depen- 
dence, a feature typical of strongly correlated systems 
H. In contrast, however, the specific heat experiments 
revealed a rather small effective mass enhancement fj). 
Recent X-ray photoemission studies on n-layer Sr doped 
systems (n=2 and oo) by Dessau et al. [0 showed ad- 
ditional puzzling aspects. For instance, in the n=oo 
compound at X = 0.4 there is a strong suppression of 
spectral weight near the Fermi energy, which correlates 
with the appearance of a pseudogap in the n=2 bi-layer 
system. Remarkably, this bi-layer system still shows a 
strong CMR effect with a rather good metallic ferromag- 
netic (FM) state at low T. 

Providing theoretical insight to understand this puz- 
zling combination is a challenging problem of strong cur- 
rent interest and the goal of the present work. One pos- 
sible approach would be to assume that the anomalous 
behaviour is due to the combined interaction of the many 
degrees of freedom which are a priori present in the per- 
ovskite manganites, such as, spin, orbital and local lattice 
distortion p,||. Here we take a different stance and by 
some reasonable simplifying assumptions, we are lead to 
a minimal model Hamiltonian that can be reliably solved. 
It will provide valuable physical insight into the nature 
of the FM metallic state. 

A rather general Hamiltonian for the low energy 
physics of the CMR manganite systems is, MM 
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where a, b = 1,2 are the orbital indexes of the e g bands, 
a is the spin, t^ b correspond to the amplitudes for nearest 

neighbours hopping, and S denotes a local spin 3/2 due 
to the three electron in the ti g orbitals The local spin 
density s of the (1 — x) conduction electrons in the e g 
bands is coupled to the latter through the Hund's rule. 
The first two terms define a ferromagnetic Kondo lat- 
tice model and give a realization of the double exchange 
mechanism (DEM) IjlOj, the third includes Jahn- Teller 
polaronic effects JJT| . These terms are widely considered 
to be responsible for much of the physics associated to the 
CMR effect that occurs at intermediate T of the order of 
the ferromagnetic transition T c . However, it was recently 
emphasized that the model also predicts phase separation 
pjfl in various regions of parameter space, which may 
play an important role for the complete understanding of 
the CMR. The last term accounts for local correlations 
due to a Hubbard-like interaction. While this interac- 
tion is recognized as an essential ingredient in every other 
transition metal oxide, it has only received very limited 
attention in manganites so far |13| |l5[ . As we shall see, 
it will play a most important role in the anomalous FM 
phase at low T. 

The exact solution of the full model ([j]) is clearly not 
possible, so we shall make a few simplifying assumptions 
which are well justified in the FM phase at low enough 
temperatures. In passing we recall that the stability of 
the FM state is due to the DEM, that favors the gain in 
kinetic energy when the system is doped with x holes per 
site |n]|l8| ■ I n a fi rs t approximation, we can neglect the 
Jahn- Teller polaronic term Hjt since the dynamical lo- 
cal distortions of the lattice, which have a notable effect 
in insulating phases and close to T c , are found to dis- 
appear when one enters the metallic FM state |19| , pfj| ,p[ . 
On the other hand, we can also safely assume the local 
spins 3/2 to be uniform and static and the magnetiza- 
tion of the (1 — x) conduction electrons to be perfectly 
aligned to them. This assumption is justified by the ex- 
perimental observation that at low enough temperatures 
(below ~200K at x ~ 0.2) the local magnetization is 
found to saturate at the classical value 101. In conse- 



quence, in (Q) we can further drop all terms containing 
operators labeled with, say, the down spin index j. At 
this point, we may take the limit of large lattice connec- 
tivity (or d — > oo) and use the dynamical mean field the- 
ory (DMFT) pTLuJ as a controlled approximation to solve 
our model. The lattice Hamiltonian is mapped onto an 
impurity problem with a self-consistency condition, and 
the resulting local action reads, 



Sice = *^2 d l(iuj n )[G 1 (iu n )] ah d h (iuj u ) 



nab 



u 



+/3— 2_j n a (iu!n)n h (-i^n) 



(2) 



na^b 



where [G (ioJ n )]a,b obeys the self-consistency condition 
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with G the local Green's function J22] For realistic e g - 
bands the off diagonal hopping elements of the matrix 
tf£° are of the same order of magnitude as the diagonal 



ones 



thus, for simplicity, we may set t 
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With 

a rotation of the impurity operators {di^cfc} — > {7 = 
(di + d 2 )/v / 2,7 / = (di - d2)/V2}, the matrix G^ x (w n ) 
becomes diagonal and the self-consistency condition now 
reads, 
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This equation tells us that in the new basis the 
7— electrons are mobile with hopping amplitude 2t, while 
the 7'— electron are localized. As the interaction term in 
(0) is rotationally invariant, the self-consistent system of 
equations (0) and (ph becomes formally identical to that 
of a Falicov-Kimball model (FKM) pfel. In this anal- 



ogy, the orbital degrees of freedom labeled by 7 and 7' 
correspond to the light and heavy spin-less electrons of 
the FKM. 

Before discussing the physics of the model in the 
present context, we would like to show that after the var- 
ious manipulations and simplifying assumptions made we 
have not lost touch with reality, and on the contrary, we 
did gain precious insight. To this end, we use density 
functional theory and obtain the realistic bandstructure 
of ferromagnetic LaSrMnOa within the local density ap- 
proximation. The bandstructure for the majority spin 
electrons is plotted in Figfil along one of the main direc- 
tions of the cubic structure. The results shown corre- 
spond to the two Mn bands with e g symmetry averaged 
over the Brillouin zone. It can be clearly observed in 
the figure that the d z 2 band is dispersive and crosses the 
Fermi energy, while the d x 2_„2 band near the Fermi en- 
ergy is essentially flat. The results of this realistic calcu- 
lation give strong support to the simplified model that we 
derived before. We have therefore gained valuable physi- 
cal insight as it will turn out that the anomalous metallic 
behavior is due to the strong scattering of fast electrons 



moving in the d z 2 channel colliding with the very slow 
electrons in the d x 2_ y 2 channel. Indeed, this is the ba- 
sic mechanism that leads to non-Fermi liquid behavior 
in the DMFT solution of the FKM ||j26j. This effect 
occurs near half filling (in the so called "pinning" region) 
and the incoherent metallic behavior can be thought of 
arising from the superposition of the upper and lower 
Hubbard bands that are split by the interaction strength 
U. 




FIG. 1. LDA bands for majority spin Mn orbitals of FM 
Lai-zA^MnOs with x=0.175. The results are for bands with 
d z 2 (a) and d x 2_ y 2 (b) character along the K z direction. The 
intensity corresponds to the average over K x and K y in the 
Brillouin zone. Doping is similuted in the Virtual Crystal 
Approximation. 

We shall now show that many anomalous properties 
of the FM metallic state observed in Sr— doped man- 
ganites are qualitatively captured by our model, using 
realistic values for the parameters. For the particular 
case of the FKM it is well known ||||6]|l|] that the 
qualitative behavior of the spectral functions does not 
depend on the choice of the non-interacting density of 
states (DOS). Therefore, we shall use a Lorentzian DOS 
P5| , [l8| that allow for fast and simple calculations. In our 
units 2t=leV, that gives a bandwidth W(FWMH)=2eV, 
as suggested by bandstructure calculations. The value of 
the local Coulomb repulsion is set to £/=3eV, which is 
consistent with the experimental estimates of either the 
charge transfer or the Mott gap in Ref. p7[ . 

In FigJ| we show the results for the optical conductiv- 
ity [g6| at various temperatures and small doping p6[ . 
We clearly observe the lack of a Drude part which is due 
to the non-Fermi liquid character of the metallic state. 
This Drude contribution is absent even at T = and we 
can understand this result as due to the presence of lo- 
calized electrons that break the translational invariance. 
More physically, the incoherent metallic state is due to 
the strong scattering of itinerant electrons at the posi- 



tions occupied by the localized ones. We thus find that 
the optical spectra is dominated by the mid-infrared con- 
tribution which has a T dependence very similar to the 
experimental results of Okimoto et al. jj| . Note that the 
lower values of T were chosen to be of the same order of 
magnitude as in the experiments. 
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FIG. 2. Optical conductivity (j{uj) calculated for various 
temperatures T/W=0.005, 0.015, 0.05, 0.15 (solid, dashed, 
long dashed, dot-dashed) and total number of particles 
n=0.93. 



metallic state looses its non-Fermi liquid character. As 
we show in the inset of the figure, the reason for this be- 
haviour is that when we reduce the chemical potential fi 
(decreasing the total number of particles) the occupation 
of the localized 7'-electrons rapidly decreases and, at the 
same time, the occupation of the itinerant 7-electrons 
increases. Thus we have the simultaneous effects of a 
reduction of the scattering and an increase in the num- 
ber of carriers with doping. We would like to point out, 
that this behavior is in good qualitative agreement with 
recent experimental results of Tokura et al. p8[ . 
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FIG. 3. Optical conductivity a(u>) for T/W = 0.005 and 
the total number of particles n =0.95 (solid), 0.93 (dashed), 
0.9 (long dashed) and 0.8 (dash-dotted). The inset shows the 
number occupation n of the 7 (solid) and 7' (dashed) orbitals 
as a function of the chemical potential /j,. 

In Fig.[| we plot the results for the optical response 
as we doped the system away from half filling. We ob- 
serve that the Drude part contribution rapidly emerges 
as the model recovers its translational invariance and the 



FIG. 4. The top figure is the density of states p(u>) ob- 
tained for T=0.005 and n=0.93. The bottom figure shows the 
e-resolved density of states p e (oj) calculated for various val- 
ues of e across the Brillouin zone. Note that the vertical axis 
correspond to the frequency u) and the the curves are plot- 
ted along the horizontal axis to better display the dispersion 
of the quasiparticle feature. The spectra p t {uj) were shifted 
(in the horizontal direction) by their corresponding e value. 
The thick dotted line is an aid to the eye to follow the dis- 
persion of the quasiparticle feature as a function of e across 
the Brillouin zone. The thin dashed line shows the itinerant 
7-electron band dispersion of the non interacting FK model. 

We now turn to the calculation of the density of states 
(DOS) which allows for the interpretation of the remark- 
able X-ray photoemission results Q. The photoemis- 
sion experiments were performed on n— layer systems at 
x = 0.4. The n~1 is a bi-layer system and the n=oo cor- 
responds to the usual 3-dimensional compound. In the 
first case the samples can be cleaved and angle resolved 
photoemission results were obtained. They revealed the 
unexpected feature of a well defined quasiparticle reso- 
nance that shows usual band dispersion behavior away 
from the Fermi surface, but the dispersion flattens out 
as it approached the (expected location of the) Fermi 
surface, and remains flat without crossing uj—0. After 
angle integration, this behavior leads to the presence 



of a so called pseudogap at the Fermi energy. In the 
n=oo compound, angle resolved experiments are not pos- 
sible and the angle integrated photoemission also shows 
a strong suppression of the spectral weight at Ep . In 
the upper part of Fig.H we show the model prediction 
for the DOS that clearly show a suppression of spectral 
weight near u>=0. This is due to the splitting of the up- 
per and lower Hubbard bands by an energy of order U. 
To gain further insight on the nature of the Hubbard 
bands and to make comparison of the predictions of the 
model to the angle resolved data, we compute the ana- 
log to a momentum resolved DOS within DMFT. In this 
theory |2^,0l the Green's function implicitly depends on 
momenta through the single particle energy e = ej: as, 
G(e, u) — l/(w — e — S(w)) where S(w) is the local self- 
energy. Therefore, the connection to the usual fc-resolved 
DOS is made by noting that momenta close to the zone 
center (r-point) correspond to values of e ~ —W; mo- 
menta near Ep, to values of e ~ 0; and momenta close to 
the zone boundary, toe~ W. In the lower part of Fig.0 
we plot the results for the imaginary part of G(e, u>) for 
various values of e across the Brillouin zone. These re- 
sults provide insights on the nature of the dispersion of 
the interacting model and can be qualitatively compared 
to the angle resolved data of Ref. Q . Using the same re- 
alistic values for the model parameters, we clearly see in 
the figure a distinct quasiparticle feature that disperses 
across the lower Hubbard band in the occupied part of 
the spectra (to < 0), and another one that has similar be- 
haviour in the upper Hubbard band for the unoccupied 
part (u) > 0). The first one, which is to be contrasted to 
photoemission experiments, shows a rather usual disper- 
sive pattern for e ~ — W (near the zone center), but it 
then flattens out as e — > 0, and remains flat up to higher 
values of e (near the zone boundary). It is important 
to note that as the quasiparticle peak fails to cross the 
Fermi energy, it also becomes broader and has a rapid 
loss of weight, indicating its progressive incoherent na- 
ture as it approaches the low frequencies u> ~ 0. This 
behavior is in remarkable qualitative agreement with the 
angle resolved data in the bi-layer system 0] . 

To conclude, we have derived and studied a model rel- 
evant for the anomalous FM metallic state of Sr doped 
manganites at low T. Our results provide a consistent 
description of several unusual observations in the exper- 
imental spectroscopy. The analogy of our model with an 
orbital Falicov-Kimball model, also explains the experi- 
mental finding that at upon doping to high values of x, 
the incoherent metal crosses-over to a regular Fermi liq- 
uid. Even more interestingly, at low values of doping, 
the analogy with the FKM suggest also an explanation 
for the driving force in the recent observation by Endoh 
et al. ]3l| of charge-ordering in the FM insulating phase 
as due to the well known charge density wave instability 
p4|,p9| of the model in bipartite lattices near half filling. 



We stress that our present work remains at a qualitative 
level, and a detailed calculation of the phase boundaries 
p0| and the eventual modifications that one may have to 
introduce to make more precise comparison with experi- 
ments are left for future work. 
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